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Abstract 

The standard model of non-relativistic quantum electrodynamics describes non-relativistic quantum 
matter, such as atoms and molecules, coupled to the quantized electromagnetic field. Within this model, we 
review basic notions, results and techniques in theory radiation. We describe the key technique in this area - 
the spectral renormalization group. Our review is based on joint works with Volker Bach and Jiirg Frohlich 
and with Walid Abou Salem, Thomas Chen, Jeremy Faupin and Marcel Griesemer. Brief discussion of 
related contributions is given at the end of these lectures. This review will appear in "Quantum Theory 
from Small to Large Scales", Lecture Notes of the Les Houches Summer Schools, volume 95, Oxford 
University Press, 201 1. 
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1 Overview 

We will describe some key results in theory of radiation for the standard model of non-relativistic elec- 
trodynamics (QED). The non-relativistic QED was proposed in early days of Quantum Mechanic^] and it 
describes quantum-mechanical particle systems coupled to quantized electromagnetic field. It arises from 
a standard quantization of the corresponding classical systems (with possible addition of internal - spin - 
degrees of freedom J] and it gives a complete and consistent account of electrons and nuclei interacting with 
electro-magnetic radiation at low energies. In fact, it accounts for all physical phenomena in QED, apart from 
vacuum polarization. Sample of issues it addresses are 

• Stability; 

• Radiation; 

• Renormalization of mass; 

• Anomalous magnetic moment; 

• One-particle states; 

• Scattering theory. 

There was a remarkable progress in the last 10 or so years in rigorous understanding of the corresponding 
phenomena. In this brief review we will deal with results concerning the first two items. We translate them 
into mathematical terms: 

*Dept. of Math., Univ. of Toronto, Toronto, Canada; Supported by NSERC Grant No. NA7901 
'it was used, as already known, by Fermi (|37|) in 1932 in his review of theory of radiation. 
2 In fact, it is the most general quantum theory obtained as a quantization of a classical system. 
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• Stability <^=> Existence of the ground state; 

• Radiation Formation of resonances out of the excited states of particle systems, scattering theory. 

One of the key notions here is that of the resonance. It gives a clear-cut mathematical description of processes 
of emission and absorption of the electro-magnetic radiation. 

The key and unifying technique we will concentrate on is the spectral renormalization group. It is easily 
combined with other techniques, e.g. complex deformations (for resonances), the Mourre estimate (for dy- 
namics), analyticity, fiber integral decompositions and Ward identities (used so far for translationally invariant 
systems). It was also extended to analysis of existence and stability of thermal states. 
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2 Non-relativistic QED 
2.1 Schrodinger equation 

We consider a system consisting of n charged particles interacting between themselves and with external 
fields, which are coupled to quantized electromagnetic field. The starting point of the non-relativistic QED 
is the state Hilbert space H — H p (£)Hf, which is the tensor product of the state spaces of the particles, H p , 
say, H p = L 2 (U. 3n ), and of Bosonic Fock space Hf of the quantized electromagnetic field, and the standard 
quantum Hamiltonian H = H gx on H = H p (£)Hf, given (in the units in which the Planck constant divided 
by 2ir and the speed of light are equal to 1 : h = 1 and c = 1) by 

n 

H = Y, ^( iV ^ ~ 9A x { Xj )f + V(x) + H f (1) 
j=i i 

(see IT371 and ||96l ). Here, rtij and Xj, j = 1, n, are the ('bare') particle masses and the particle positions, 
x = (xi, . . . ,!„), V(x) is the total potential affecting particles and g > is a coupling constant related 
to the particle charges, A x :— \ * A, where A(y) is the quantized vector potential, in the Coulomb gauge 
(drv A(y) = 0), describing the quantized electromagnetic field, and \ is an ultraviolet cut-off, 

A X (V) = / (e iky a(k) + e-* fe %*(A0) X (fc)^= (2) 

J y/\k\ 



(a(k) and a*(k) are annihilation and creation operators acting on the Fock space Hf = T , see Supplement 
[Hjfor the definitions), and Hf is the quantum Hamiltonian of the quantized electromagnetic field, describing 
the dynamics of the latter, it is given by 

Hj = J d 3 ku}(k)a*(k)-a(k), (3) 

where w(fc) = \k\ is the dispersion law connecting the energy of the field quantum with its wave vector k. 
For simplicity we omitted the interaction of the spin with magnetic field. (For a discussion of this Hamiltonian 
including units, the removal of the center-of-mass motion of the particle system and taking into account the 
spin of the particles, see Appendix [A] Note that our units are not dimensionless. We use this units since we 
want to keep track of the particle masses. To pass to the dimensionless units we would have to set m e \ = 1 
also.) The Hamiltonian H determines the dynamics via the time-dependent Schrodinger equation 



idtip = Hip, 
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where tp is a differentiable path in H = "H p ®Tif. 

The ultraviolet cut-off, x, satisfies x(A:) = 1 in a neighborhood of k = and is decaying at infinity on 
the scale k and sufficiently fast. We assume that V(x) is a generalized n-body potential, i.e. it satisfies the 
assumptions: 

(V) V(x) = J2i Wi(iTix), where 7Tj are a linear maps from M 3 ™ to M. mi , m.j < 3n and Wi are Kato-Rellich 
potentials (i.e. W^x) € LP<(ffi m «) + (L° c (R 3n )) e withp, = 2 for m, < 3, p 4 > 2 for to, = 4 and 
Pi > rrii/2 for > 4). 

Under the assumption (V), the operator H is self-adjoint and bounded below. 

We assume for simplicity that our matter consists of electrons and the nuclei and that the nuclei are 
infinitely heavy and therefore are manifested through the interactions only (put differently, the molecules are 
treated in the Born - Oppenheimer approximation). In this case, the coupling constant g is related to the 
electron charge — e as g :— o?l 2 , where a = « j^, the fine-structure constant, and mj = m. It is 
shown (see Section [T2| and a review in Q) for references and discussion) that the physical electron mass, 
rrieh is not the same as the parameter rn = mj (the 'bare' electron mass) entering ([TJ, but depends on m and 
k. Inverting this relation, we can think of m as a function of m c \ and k. If we fix the particle potential V{x) 
(e.g. taking it to be the total Coulomb potential), and m e \ and e, then the Hamiltonian ([1} depends on one 
free parameter, the bare electron mass m (or the ultraviolet cut-off scale, k). 



2.2 Stability and radiation 

We begin with considering the matter system alone. As was mentioned above, its state space, T~L p , is ei- 
ther L 2 (M. 3n ) or a subspace of this space determined by a symmetry group of the particle system, and its 
Hamiltonian operator, H p , acting on W p , is given by 



where is the Laplacian in the variable Xj and, recall, V(x) is the total potential of the particle system. 
Under the conditions on the potentials V(x), given above, the operator H p is self-adjoint and bounded below. 
Typically, according to the HVZ theorem, its spectrum consists of isolated eigenvalues, e < e{ < ... < 
£( p ), and continuum [E^ p ', oo), starting at the ionization threshold £^0, as shown in the figure below. 



The eigenfunctions corresponding to the isolated eigenvalues are exponentially localized. Thus left on its 
own the particle system, either in its ground state or in one of the excited states, is stable and well localized 
in space. We expect that this picture changes dramatically when the total system (the universe) also includes 
the electromagnetic field, which at this level must be considered to be quantum. As was already indicated 
above what we expect is the following 

• The stability of the system under consideration is equivalent to the statement of existence of the ground 
state of H, i.e. an eigenfunction with the smallest possible energy. 




(4) 



Hidden Picture 



• The physical phenomenon of radiation is expressed mathematically as emergence of resonances out of 
excited states of a particle system due to coupling of this system to the quantum electro-magnetic field. 
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Our goal is to develop the spectral theory of the Hamiltonian H and relate to the properties of the relevant 
evolution. Namely, we would like to show that 

1) The ground state of the particle system is stable when the coupling is turned on, while 

2) The excited states, genetically, are not. They turn into resonances. 



2.3 Ultra-violet cut-off 

We reintroduce the Planck constant, h, speed of light, c, and electron mass, m c \, for a moment. Assuming the 
ultra-violet cut-off decays on the scale k, in order to correctly describe the phenomena of interest, such 
as emission and absorption of electromagnetic radiation, i.e. for optical and rf modes, we have to assume that 
the cut-off energy, 

hen ^> a 2 m c \c 2 , ionization energy, characteristic energy of the particle motion. 

On the other hand, we should exclude the energies where the relativistic effects, such as electron-positron 
pair creation, vacuum polarization and relativistic recoil, take place, and therefore we assume 

hen -C to c ic 2 , the rest energy of the the electron. 

Combining the last two conditions we arrive at a 2 m c \c/h <C k <C m e \c/K, or in our units, 

a 2 m e i <C n -C m e i . 

The Hamiltonian ([T} is obtained by the rescaling x — > a x and k — > a 2 k of the original QED Hamiltonian 
(see Appendix [a|. After this rescaling, the new cut-off momentum scale, k' — a~ 2 K, satisfies 

m c \ <C k' <C a~ 2 m e \, 

which is easily accommodated by our estimates (e.g. we can have k' = 0(a~ 1 ^ 3 m e i)). 



3 Resonances 

As was mentioned above, the mathematical language which describes the physical phenomenon of radiation 
is that of quantum resonances. We expect that the latter emerge out of excited states of a particle system due 
to coupling of this system to the quantum electro-magnetic field. 
Quantum resonances manifest themselves in three different ways: 

1) Eigenvalues of complexly deformed Hamiltonian; 

2) Poles of the meromorphic continuation of the resolvent across the continuous spectrum; 

3) Metastable states. 

3.1 Complex deformation 

To define resonances we use complex deformation method. In order to be able to apply this method we 
choose the ultraviolet cut-off, x{k), so that 

The function 9 —> x( e ~ 9 k) has an analytic continuation from the real axis, K, to the strip {0 £ 
C||Im 0\ < tt/4} as a L 2 f]L°°(R 3 ) function, 
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e.g. x(k) = e ' fe ' l K . For the same purpose, we assume that the potential, V(x), satisfies the condition: 

(DA) The the particle potential V(x) is dilation analytic in the sense that the operator-function 9 — > V(e x) 
(—A + has an analytic continuation from the real axis, R, to the strip {9 € C||Im 9\ < 9o} for 
some 9 n > 0. 

To define the resonances for the Hamiltonian H we pass to the one-parameter (deformation) family 

H e := UgHUg 1 , (5) 

where 9 is a real parameter and Ug, on the total Hilbert space T~L := W p ® T, is the one-parameter group of 
unitary operators, whose action is rescaling particle positions and of photon momenta: 

Xj — > e 8 Xj and k — > e~ e fc. 

One can show show that: 

1) Under a certain analyticity condition on coupling functions, the family Hg has an analytic continuation 
in 9 to the disc D(0, 9 ), as a type A family in the sense of Kato; 

2) The real eigenvalues of Hg, lm9 > 0, coincide with eigenvalues of H and that complex eigenvalues 
of Hg, lm9 > 0, lie in the complex half-plane C~; 

3) The complex eigenvalues of Hg, lm9 > 0, are locally independent of 9. The typical spectrum of 
[Hg = Hg]\ g= o, lm9 > (here the superindex SM stands for the standard model) is shown in the figure 
below. 

Hidden Picture 



We call complex eigenvalues of Hg, Im 9 > the resonances of H. 

As an example of the above procedure we consider the complex deformation of the hydrogen atom and 
photon Hamiltonians H^y^ ;= — — At and Hf. 

tj -20/ ^ \\ -8 a tj -Sit 

H hydr g = e (- — A)-e j-r, H fe = e H f . 

Let e^ ydr be the eigenvalues of the hydrogen atom. Then the spectra of these deformations are 

a{H hydr0 ) = {e'; ydr } U e- 2Im >,oo), a(H f6 ) = {0} U e - Ime [0,oo). 



3.2 Resonances as poles 

Similarly to eigenvalues, we would like to characterize the resonances in terms of poles of matrix elements 
of the resolvent (H — z)^ 1 of the Hamiltonian H. To this end we have to go beyond the spectral analysis of 
H. Let ^g = Ug^, etc., for 9 e K and z e C+. Use the unitarity of Ug for real 9, to obtain (the Combes 
argument) 

(*, (H-z)- 1 $) = {y g ,(Hg-z)- 1 $ e )- (6) 

Assume now that for a dense set of 'f's and $'s (say, T>, defined below), ^g and $e have analytic continua- 
tions into a complex neighbourhood of 9 = and continue the r.h.s of |6]) analytically first in 9 into the upper 
half-plane and then in z across the continuous spectrum. This meromorphic continuation has the following 
properties: 

• The real eigenvalues of Hg give real poles of the r.h.s. of |6| and therefore they are the eigenvalues of 
H. 
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• The complex eigenvalues of Hg are poles of the meromorphic continuation of the l.h.s. of |6]) across 
the spectrum of H onto the second Riemann sheet. 

The poles manifest themselves physically as bumps in the scattering cross-section or poles in the scattering 
matrix. 

The r.h.s. of (|6]l has an analytic continuation into a complex neighbourhood of 9 = 0, if $ £ V, where 

V := (J Ra,n(xN<nX\T\<o)- (7) 

n>0 : a>0 

Here N — J d 3 ka* (k)a(k) be the photon number operator and T be the self-adjoint generator of the one- 
parameter group Ug, (It is dense, since N and T commute.) 



3.3 Resonance states as metastable states 

While bound states are stationary solutions, one expects that resonances to lead to almost stationary, long- 
living solutions. Let z*, Imz, < 0, be the ground state or resonance eigenvalue. One expects that for an 
initial condition, ipo, localized in a small energy interval around the ground state or resonance energy, Re z*, 
the solution, ip — e~ lHt ipo, of the time-dependent Schrodinger equation, idtip = Hip, is of the form 

^ = e -«*V* + loc (t- Q ) + O res (/), (8) 

for some a, (3 > (depending on ipo), where 

• (j)* is either the ground state or an excited state of the unperturbed system, depending on whether z* is 
the ground state energy or a resonance eigenvalue; 

• The error term Oi oc (^ Q ) satisfies ||(1 + \T\)- u Oi oc (t- a )\\ < Ct- a , where T is the generator of the 
group Ug, with an appropriate v > 0. 

For the ground state, ([8j, without the error term O les (g^), is called the local decay property (see Section 
10 1. One way to prove it is to use the formula connecting the propagator and the resolvent: 

1 f°° 

f(H) = - / dXf(X)e- lXt lm (if — A — iQ)- 1 . (9) 

7T J-oo 



e -mt 



Then one controls the boundary values of the resolvent on the spectrum (the corresponding result is called 



the limiting absorption principle, see Appendix D.2 1 and uses properties of the Fourier transform 



z 



For the resonances, (|8]l implies that — Im has the meaning of the decay probability per unit time, and 
(-Imz,) -1 , as the life-time of the resonance. To prove it, one uses |6]) and the analyticity of its r.h.s. in 
and performs in Q a suitable deformation of the contour of integration to the second Riemann sheet to pick 
up the contribution of poles there. This works when the resonances are isolated. In the present case, they are 
not. This is a consequence of the infrared problem. Hence, determining the long-time behaviour of e~ lHt ipo 
is a subtle problem in this case. 



3.4 Comparison with Quantum Mechanics 

This situation is quite different from the one in Quantum Mechanics (e.g. Stark effect or tunneling decay) 
where the resonances are isolated eigenvalues of complexly deformed Hamiltonians. This makes the proof 
of their existence and establishing their properties, e.g. independence of 9 (and, in fact, of the transformation 
group Ug), relatively easy. In the non-relativistic QED (and other massless theories), giving meaning of the 
resonance poles and proving independence of their location of 9 is a rather involved matter. 
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3.5 Infrared problem 

The resonances arise from the eigenvalues of the non-interacting Hamiltonian H g=0 . The latter is of the form 

H = Hpart ® 1/ + Ipart ® #f ■ (10) 

The low energy spectrum of the operator H consists of branches [e| p ' , oo) of absolutely continuous spectrum 
and of the eigenvalues e^'s, sitting at the continuous spectrum 'thresholds' e| p ''s. Here, recall, < < 

... < £( p ) are the isolated eigenvalues of the particle Hamiltonian H p . Let be the eigenfunctions of the 

(v) 

particle system, while be the photon vacuum. The eigenvalues e\ • s correspond to the eigenfunctions 

<j)f^ (g> Q of TJo- The branches [e\ p \ oo) of absolutely continuous spectrum are associated with generalized 

(v) 

eigenfunctions of the form <p\ ' ®g\, where g\ are the generalized eigenfunctions of Hf : Hfg\ = Xg\, < 
A < oo. 

The absence of gaps between the eigenvalues and thresholds is a consequence of the fact that the photons 
are massless. To address this problem we use the spectral renormalization group (RG). The problem here is 
that the leading part of the perturbation in H is marginal. 



4 Existence of the Ground and Resonance States 
4.1 Bifurcation of eigenvalues and resonances 

Stated informally what we show is 

• The ground state of H\ g=0 the ground state H (e = + 0(g 2 ) and e < 4 P) ); 

• The excited states of H\ g=0 =>■ (generically) the resonances of H (ej t k = + 0(g 2 )); 

• There is E > inf <j(H) (the ionization threshold, E + E( p ) + 0(g 2 )) s.t. for energies < E that particles 
are exponentially localized around the common center of mass. 

For energies > E the system either sheds off locally the excess of energy and descends into a localized state 
or breaks apart with some of the particles flying off to infinity. 
To formulate this result more precisely, denote 

e g %(u) := min{|eW - e« | | i * j, e^,ef < „}. 

Theorem 4.1 (Fate of particle bound states). Fix e^f' < v < inf a ess (H p ) and let g <C £gap{v). Then 
for 5 ^0, 

• H has a ground state, originating from a ground state of H\ g=0 (e = + 0(g 2 ), e < e^); 

• Generically, H has no other bound state (besides the ground state); 

• Eigenvalues, e^ p ' < u, j ; ^ 0, of H\ g=0 => resonance eigenvalues, Cj t k, of H; 

• e j,k = ^ + 0(g 2 ) and the total multiplicity of ej^ equals the multiplicity of e^; 

• The ground and resonance states are exponentially localized in the physical space: 

\\e 5 ^7p\\ < oo, VV> G R&nE A (H), S < - sup A. 
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Hidden Picture 

(v) 

Remark. The relation e < e o i s due to the fact that the electron surrounded by clouds of photons become 
heavier. 

4.2 Meromorphic continuation across spectrum 

Theorem 4.2. (Meromorphic continuation of the matrix elements of the resolvent) Assume g -C ef} p {i>) and 
let eo := inf a(H) be the ground state energy of H. Then 

• For a dense set (defined in |7]) below) of vectors ^ and $, the matrix elements 

F(z,*,$) := (H - z)- 1 ®) 
have meromorphic continuations from C + across the interval (eo, v) C a ess (H) into 

{zeC-\e <Rez <is}/ (J S jtk , 

0<j<j(u) 

where Sj k are the wedges starting at the resonances 

S jtk := {z e C | ^Re (e e (z - e j>k )) > |Im (e e (z - e i>fe ))|}; (11) 

• This continuation has poles at e^: lim z _j. ej fc (e^.t — z)F(z, ^E 1 , $) is finite and ^ 0. 



4.3 Discussion 

• Generically, excited states turn into the resonances, not bound states. 

• The second theorem implies the absolute continuity of the spectrum and its proof gives also the limiting 
absorption principle for H (see Appendix D.2 for the definitions). 

• The proof of first theorem gives fast convergent expressions in the coupling constant g for the ground 
state energy and resonances. 



• One can show analyticity of tj k in the coupling constant g (see Appendix 10 for a result on the ground 
state energies). 

• The meromorphic continuation in question is constructed in terms of matrix elements of the resolvent 
of a complex deformation, He-, Im# > 0, of the Hamiltonian H. 



• A description of resonance poles is given in Section 10 
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4.4 Approach 

The main steps in our analysis of the spectral structure of the quantum Hamiltonian H are: 

• Perform a new canonical transformation (a generalized Pauli-Fierz transform) 

H — > H PF :=e- i9F He igF , 
in order to bring H to a more convenient for our analysis form; 

• Apply the spectral renormalization group (RG) on new - momentum anisotropic - Banach spaces. 

The main ideas of the spectral RG are as follows: 

• Pass from a single operator H PF to a Banach space B of Hamiltonian-type operators; 

• Construct a map, 1Z P , (RG transformation) on B, with the following properties: 

(a) 1Z P is 'isospectral'; 

(b) 1Z P removes the photon degrees of freedom related to energies > p. 

• Relate the dynamics of semi-flow, IZ™, n > 1, (called renormalization group) to spectral properties of 
individual operators in B. 



5 Generalized Pauli-Fierz transformation 

We perform a canonical transformation (generalized Pauli-Fierz transform) of H in order to bring it to a form 
which is accessible to spectral renormalization group (it removes the marginal operators). For simplicity, 
consider one particle of mass 1. We define the generalized Pauli-Fierz transformation as: 

h pf ._ e -igF He tgF^ (12) 
where F(x) is the self-adjoint operator given by 

f(z) = E /(/,,A(fcK(fc) + /,, A (fcK(fc))^p, (13) 

with the coupling function f x ,\{k) chosen as 

/^W---*!^, (14) 

v\ k \ 



with ip e C 2 bounded, with bounded derivatives and satisfying (p'(0) — 1. For the standard Pauli-Fierz 
transformation, we have ip(s) = s. 

The Hamiltonian H PF is of the same form as H. Indeed, using the commutator expansion 

e -igF(x) Hfe igF(x) = _ g 2 ^ ^ jj^ wg compute 

H PF = 7^(P + gA xv (x)) 2 + V g (x) +H f + gG(x), (15) 
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where 

A x(p (x) = ^2 J {<Px,\(k)a\(k) + tp Xt \(k)a* x (k))* ' 



A 



with the new coupling function ip\, x (k) := e\(k)e tkx — V x f x .\(k) and 

V g (x) :=V(x) + 2g 2 J2 [ \k\\f x ,x(k)\ 2 d 3 k, 

A J 

G(x):=-iJ2 I \k\(f x ,x(k)a x (k) - f x ^k)a\{k))^^. 



A " 

The potential V g (x) is a small perturbation of V(x) and the operator G(x) is easy to control. The new 
coupling function has better infrared behaviour for bounded \x\; 

\fx,x{k)\ < const min(l, <\/\k\{x)). (16) 

To prove the results above we first establish the spectral properties of the generalized Pauli-Fierz Hamil- 
tonian H PF and then transfer the obtained information to the original Hamiltonian H. 



6 Renormalization Group Map 

The renormalization map is defined on Hamiltonians acting on T~L / which as follows 

K p = p-'Sp o F p , (17) 
where p > 0, S p : B[H] — > B[H] is the scaling transformation: 

S p (l) := 1, S p (a*(k)):= p-V*a*(p- l k), (18) 
and F p is the Feshbach-Schur map, or decimation, map, 

F p {H) := Xp (H - Hx p {x p Hx p )- l X p H)x P , (19) 
where \ p and \ p is a pair of orthogonal projections, defined as 

X P = XH v9 =e 3 ® XH f < P and x P 1 - X P - 

Remark. For simplicity we defined the decimation map as the Feshbach-Schur map. Technically it is more 
convenient to use the smooth Feshbach-Schur map, which we defined and discussed in Appendix [D] The 
smooth Feshbach-Schur map uses 'smooth' projections which form a partition of unity \ 2 p + x% — 1> instead 
of true projections as defined above. 

The map F p is isospectral in the sense of the following theorem: 
Theorem 6.1. (i) A € p(H) «0e p(F p (H - A)); 

(ii) H^ = \ii ^ F P {H - A) tp = 0; 

(iii) dimNull(iJ - A) = dimNullFp(H - A); 

(iv) (H -X)- 1 exists <^ F P (H — A) -1 exists. 

For the proof of this theorem as well as for the relation between ip and ip in (ii) and between (H — A) -1 
and F p (H — A) -1 in (iv) see Appendix [c] 
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7 A Banach Space of Hamiltonians 

We will study operators on the subspace Ranxi of the Fock space T . Such operators are said to be in the 
generalized normal form if they can be written as: 

m+n>0 

/m+n m m-\-n 

J\ d 3 ki Y\_a*(ki) W mi n(Hf,k( m+n )) Y[ a ( k i): 

3 T " i=l j=l i=m+l 

where B\ denotes the Cartesian product of r unit balls in M 3 , k( m ) '■= (ki, . ■ ■ , k m ) and w m , n : I x 
j^m+n _^ ^ j ._ jo, 1]. We sometimes we display the dependence of H and W,„ „ on the coupling 
functions w := (u> m , n , m + n > 0) by writing H[w] and W m ,„[W]. 

We assume that the functions ti! m , ra (r, , k( m + n \) are continuously differentiable in r € 7, symmetric 

w.r.t.thevariables(fci,...,fc m )and(fc m+ i,...,fc m+n )andobey ||w m ,„|| M ,i := Z)* =0 ||9?«WilU < °o, 
where /x > and 

m+n 

||«VnlU max sup \\k \~^ J\ \ k i\ 1/2w ™,n(r; k( m+n ))\. (21) 

3 rS7,fe ( „ l+ „ ) eB 1 m+ " i=1 

Here fej g M 3 is the j— th 3— vector in fe( m>n ) over which we take the supremum. Note that these norms are 
anisotropic in the total momentum space. 

For /i > and < £ < 1 we define the Banach space 

&*:={H : ||fr|| Mi£ := £ r (m+n) ||«V«ILi < oo}. (22) 

m+n>0 

We mention some properties of these spaces 

• For any p > and < £ < 1, the map H : w—> H[w], given in ( [20] ), is one-to-one. 

• If If is self-adjoint, then so are W ,o and J2 m +n>i XiW m ,nXi ( see p0|>)- 

Remarks. 1) Unlike the Banach spaces defined in iTPJl [141 [7J, the Banach spaces are anisotropic in the 
momentum space. This is needed to overcome the problem of marginal operators which arise in the renor- 
malization group approach. 

2) The self-adjointness statement follows from [7], Eq. (3.33). 



7.1 Basic bound 

The following bound shows that our Banach space norm control the operator norm and the terms with higher 
numbers of creation and annihilation operators make progressively smaller contributions: 

Theorem 7.1. Let Xp = Xh s < p - Then for all p > and m + n> 1 

pin+n+p 

\\XpHm,nXp\\ < — 7=^=T II W m,n II ft- (23) 
V to! n\ 

Sketch of proof. For simplicity we prove this inequality for m = n = 1. Let <f> E T and $j. = a(k)x P 4>- We 
have 

r 2 

( Xp <i>,H ltl Xp<P) = / Y[d 3 h ($ kll w 1>1 (Hf,k 1 ,k 2 )®k 2 }- 
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Now we write \p — XpX2p an d pull Xip toward w\ \{Hf\ fci, k 2 ) using the pull-trough formulae 

a(k)f(H f ) = f(H f + \k\)a(k), f(H f )a*(k) = a* (k) f (H f + \k\) 
(see Appendix [G]>. This gives 

\{<t>,Xp H i,ix P 4>)\ 

2 

FT d 3 fc, ; ($ fcl , x 2 p- 1 fci I wi.i (iJ/ ; fci , fc 2 ) X2 P - |fc 2 1 $fe 2 ) i 

|fei|<2p,i=l,2 i=1 



< 



< 



IJd 3 ** ||* fcl || |Ki(ff/;fci,* a )|| Pfe 2 

|fci|<2p,»=l,2 j =1 



|fci|<2/3,i=l,2 i=1 



l*ll|*2| 



Now, using |K! (H f ;k 1 ,k 2 )\\ < \\w hl \\ ^7 lk ^ and J d 3 k\\ ^/\k\<f> k \\ 2 

\ki\ 2 \k 2 \ 2 

|(0,xp^i,iXp^)I <p 2+ 1KiIL- 



d 3 fc||v^$ fc || 2 . 

^/_ff/Xp0|| 2 , we find 



7.2 Unstable, neutral and stable components 

We decompose H e B^ into the components E := (£l,HSl), T := H 0fi -{Cl,H£l), W := E m +„>i 
so that 

H = El + T + W. (24) 

If we assume sup re r 0iOO \ |T"(r) — 1| <C 1, then we have T ~ Hf. These Hamiltonian components scale as 
follows 

• p _1 5 p (-ff/) = Hf (Hf is a fixed point of p~ 1 S p ); 

• p~ 1 S p (E ■ 1) = p~ x E ■ 1 (E ■ 1 expand under p~ 1 S p at a rate 

• ||5p(W min )|| M < p Q ||w m ,n|| M , a := m + n - 1 + /i# m +„=i (W mn contract under p~ 1 S p , if p > 0). 

Thus for > 0, E, T, W behave, in the terminology of the renormalization group approach, as relevant, 
marginal, and irrelevant operators, respectively. For p = 0, the operators W mn , rn + n = 1, become 
marginal. 



8 Action of Renormalization Map 

To control the components E, T, W of H we introduce, for a, 0, 7 > 0, the following polydisc: 

£>"(a,j9,7) := [i? = E + T + W e \ \E\ < a, 
sup |T'(r) - 1| < A little < 7}- 

re[0,oo) - 1 

(Strictly speaking we should write XH pg =e j ® T> ,J, (a, /3, 7) instead of T> fl (a, f3, 7) (for various sets of param- 
eters a, j3, 7) in the statement below.) 
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Theorem 8.1. Let < p < 1/2, a, /3, 7 < p/8 and p = 1/2. Then there is c > 0, s.t. 

• K P (HP F ) e P"°(ao, A),7o), a = c 5 V 0-2 , A) = c 5 V 0_1 . 7o = cgp%, 
provided g <C 1; 

• 2?^ (a, /3, 7) C D(K P ), provided p > 0; 

• "ftp : 2? M (a, /3, 7) -» C(a', /?', 7'), continuously, with 

a' = p^a + c ( 7 2 /2p) ,fi = p + c {l 2 /2p) , 7' = cp" 7 . 

Sketch of proof of the second and third properties. 1) 2? M (p/8, 1/8, p/8) C D(lZ p ). Since T'F := H — E — T 
defines a bounded operator on J 7 , we only need to check the invertibility of H TXp on Ranx p - The operator 
E + T is invertible on Ran% p : for all r € [3p/4, 00) 

ReT(r)+ReE > r - \T(r) - r\ - \E\ 

> r(l - sup|T'(r) - 1|) - \E\ 

r 

> ^(1-1/8) - e. > p 

- 4 v ' ' 8 - 2 

=> E + T is invertible and \\(E + T)~ 1 \\ < 2/ p. 

Now, by the basic estimate, |[W|| < p/8 and therefore, 

Wx.WXpiE + T)- 1 ]] <l/4 
=>■ + T + x p Wx p is invertible on Ran x p 
=> V»(p/8,1/&,P/8)CD(F P ) = D(K P ). 

2) n p : P^(a,/3,7) -> 2?^ (a', /?', 7') (normal form of TZ p (H)). Recall that x P = X-H><p and x p := 
1 - x P - Let i7 := -E + T, so that H = H Q + W. We have shown above 

||#0 -1 XJ < ^ and ||VK|| < |. 

In the Feshbach-Schur map, F p , 

F p (H) = Xp(Ho + W~WXp(Xp(Ho + W)x p y 1 XpW)xp, 
we expand the resolvent (x p (Hq + W)x p ) 1 m the norm convergent Neumann series 

oo 

F p {H) = x P [H a + Y J i- l ) SW { H ^X 2 pW) s ]x P . 

s=0 

Next, we transform the right side to the generalized normal form using generalized Wick's theorem. 

Generalized Wick's theorem. To write the product Wi^H^x^ W) in the generalized normal form 
we pull the annihilation operators, a, to the right and the creation operators, a*, to the left, apart from those 
which enter Hj. We use the rules (see Appendix |G]i: 

a(k)a*(k') = a*(k')a{k) + 5{k - k'), 

a(k)f(H f ) = f(H f + \k\)a(k), f{H f )a*{k) = a* (k) f (H f + \k\) . 
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Some of the creation and annihilation operators reach the extreme left and right positions, while the remaining 
ones contract (see the figure below). The terms with m creation operators on the left and n annihilation 
operators on the right contribute to the (m, n)— formfactor, Wm,n, of the operator W (Hq 1 ^ 2 , W} s . As the 
result we obtain the generalized normal form of F p (H): 

F p (H) = Y, W 'm,n- 
m+n>0 

The term W' 0fi = (W$)n + (W$ - (W$)n) contributes the corrections to E + T. 

Hidden Picture 



This is the standard way for proving the Wick theorem, taking into account the presence of Hf— depen- 
dent factors. (See lfT3ll for a different, more formal proof.) 

Estimating formfactors. The problem here is that the number of terms generated by various contractions 
is O(sl). Therefore a simple majoration of the series for the (to, n)— formfactor, Wm] n , of the operator 
W{Hq X x 2 p W) s will diverge badly. 

To overcome this we re-sum the series by, roughly, representing the sum over all contractions, for a given 
to and n, as 

w^ n ~(n,{w(H^x 2 P w ) s Unn), 

where \W [Hq 1 ^ WY] m ^ n is W [Hq 1 ^, W) s , with m escaping creation operators and n escaping anni- 
hilation operators deleted. Now the estimate of Wm'n is straightforward and can be written (symbolically) 

as 

£ \\x P > w'xA s+ \ 

( for the operator norm, and similarly for B^^— norm. 



9 Renormalization Group 

To analyze spectral properties of individual operators in we use the discrete semi-flow, 7£", n > 1 (called 
renormalization group) , generated by the renormalization transformation, 1Z P . By Theorem 8.1 in order to 
iterate 1Z P we have to control the expanding direction: 7Z P ((1) = /? _1 C1- To control this direction, we adjust, 
inductively, at each step the constant component (H)q :— (57, HQ) of the initial Hamiltonian, H: 

|(i?)n-e„_ 1 |<i /? "+ 1 , 



so that 



-i is the unique zero of the function A -> (7^" _1 (i? - (H) n + A)) s 
H e the domain of 



p' 

This way one adjusts the initial conditions closer and closer to the stable manifold Ai s — r\ n D(lZ 7 p l ) 

Hidden Picture 



The procedure above leads to the following results: 
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• TVi has the fixed-point manifold M. f p := CHf\ 

• TvL™ has an unstable manifold A4 U := CI; 

• 7?.™ has a (complex) co-dimension 1 stable manifold A4 S for Ai f p ; 

• M. s is foliated by (complex) co-dimension 2 stable manifolds for each fixed point. 

Hidden Picture 



Define the polydisc T> s := P M (0, /3q, 7o)> with [i > 0, (3q ~ cg 2 p^" 1 and 70 = cgp^, the same as in 



Theorem 8.1 Using the above results, we draw the following conclusions about the spectrum of an operator 
H e V s . Find A e a neighbourhood of 0, s.t. H(X) := H-X e T>{R%) and define (A) := W l p {H{\)). 
We illustrate this as putting i? (A) through the RG (black) box: 

H(\)=> RG BOX =*> H {n) (X). 

Use that for n sufficiently large, H^ n ' (A) w C^/' f° r some ( G C, Re £ > 0, to find wanted spectral infor- 
mation about H( n ' (A) in a neighbourhood of 0. Then, use the isospectrality of 7Z P to pass this information 
up the chain. These steps are described schematically as follows: 

Spectral information about H^(X) (by 'isospectrality' of TZ p ) 
Spectral information about i?(" -1 )(A) 

=> Spectral information about H(X). 
We sum up this as 

Spec info (H) <^= RG BOX <s= Spec info (H {n \X)). 

If we start with the Hamiltonian Hg, then we are interested in, we derive this way the required spectral 
information about it. 



10 Related Results 

Existence of the ionization threshold. There is E = E^ p ' + 0(g 2 ) > inf a(H) (the ionization threshold) 
s.t. for any energy interval in A C (inf <r(H), E), 

He^H < 00, W e Ran£ A (#), 8 < E - sup A. 



Analyticity. If the ground state energy e g of H is non-degenerate, then it is analytic in g and has the following 
expansion 

00 
3=1 

where are smooth function of a (recall that g :— a 3 / 2 ). 

The proof relies on the renormalization group analysis together with the analyticity and the form-factor 
rotation symmetry transfer by the Feshbach-Schur map (see Appendix |D.3[ ) and on treating two sources of 
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the dependence of H on the coupling constant g - in the prefactor of A x (y) and in its argument (see the line 
after pT)), differently, i.e. by considering the following family of Hamiltonians 




and consider g and a as independent variables. The powers of g in ( p5| come from g in ( |26] l. 

Analyticity of all parts of H. Suppose that A H> H\ = H{w x ) is of the form |20) and is analytic in 
A € S C C and that H(w x ) belongs to some polydisc T>(a,f3,j) for all A <E S. Then A H> := 
W ^ o (0), T A := < (i?/) - <o(0), W A := H x -E x - T x are analytic in A G 5. 

Resonance poles. Can we make sense of the resonance poles in the present context? Let 

Q:={z£C-\e <Rez<v}/ [j S hk . 

Theorem 10.1. For each "J and $ from a dense set of vectors (say, |7]i), the meromorphic continuation, 
F(z, ^, $), of the matrix element (H — is of the following for ni near the resonance €a of H\ 

F(z, $) = (e ilfc - z)" 1 ^, $) + r(z, $). (27) 

Here p and r(z) are sesquilinear forms in ^ and $, s.t. 

• r(z) is analytic in Q and bounded on the intersection of a neighbourhood of ej,k with Q as 

|r(z, *,$)| < Cq,^\e j:k - z|~ 7 for some 7 < 1; 

• p 7^ at least for one pair of vectors '5 and $ and p = for a dense set of vectors ^ and <f> in a finite 
co-dimension subspace. 

Local decay. For any compactly supported function /(A) with supp/ C (inf H, oo)/(a neighbourhood of 
S), and for 8 > |, i/ < ^ — g, we have that 

\\(Y)- e e - iHt f(H)(Y)- e \\ <Cr u . (28) 
Here (Y) := (1 + Y 2 ) 1 / 2 , and y denotes the self-adjoint operator on Fock space T of photon co-ordinate, 

Y := J d 3 k a*(k)iV k a(k), (29) 

extended to the Hilbert space % = "H p ® J 7 . (A self-adjoint operator H obeying ( |28] i is said to have the 
(A, v 1 Y, 0) - local decay (LD) property.) ( |28| i shows that for well-prepared initial conditions i/jq, the proba- 
bility of finding photons within a ball of an arbitrary radius R < oo, centered say at the center-of-mass of the 
particle system, tends to 0, as time t tends to oo. 

Note that one can also show that as t — > oo, the photon coordinate and wave vectors in the support of the 
solution e~ lHt ipo of the Schrodinger equation become more and more parallel. This follows from the local 
decay for the self-adjoint generator of dilatations on Fock space T, 

B := ~ J d 3 ka*{k) {/c-Vfc+Vfc-fc} a(k). (30) 

(In fact, one first proves the local decay property for B and then transfers it to the photon co-ordinate operator 
Y.) 
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11 Conclusion 

Apart from the vacuum polarization, the non-relativistic QED provides a good qualitative description of the 
physical phenomena related to the interaction of quantized electrons and nuclei and the electro-magnetic field. 
(Though construction of the scattering theory is not yet completed and the correction to the gyromagnetic ratio 
is not established, it is fair to conjecture that while both are difficult problems, they should go through without 
a hitch.) 

The quantitative results though are still missing. Does the free parameter, m (or k), suffice to give a 
good fit with the experimental data say on the radiative corrections? Another important open question is the 
behaviour of the theory in the ultra-violet cut off. 

12 Comments on Literature 

These lectures follow the papers II 1031 1451 [TI which in turn extend HUE] El. The papers lfT03ll45l [Tl use the 
smooth Feshbah-Schur map (0[3T]), which is much more powerful (see Appendix [D|, while in these lectures 
we use, for simplicity, the the original, Feshbach-Schur map (see |fT3l[T4l ). which is simpler to formulate. 

The self-adjointness of H is not difficult and was proven in 031 for sufficiently small coupling constant 
(g) and in f72l (see also 1701 ). for an arbitrary one. 

Theorems 4.1 and 4.2 were proven in lfT3lfT4l[T5ll for 'confined particles' (the exact conditions are some- 
what technical) and in the present form in II 1031 . 

The results of |fT31l on existence of the ground state were considerably improved in [68 69l ITTl 1771 |4l |89l 
(by compactness techniques) and [9| (by multiscale techniques), with the sharpest result given in |[53~1 [87l . 
(The papers |[T5l [68] |53l [85] [86] |87l [771 HU include the interaction of the spin with magnetic field in the 
Hamiltonian.) 

Related results: 

The asymptotic stability of the ground state (local decay, see SectionfTO}: fTTl l44l l46l . 
The survival probabilities of excited states (see |8]l): |[T5l l62l [T1 . 
Atoms with dynamic nuclei: 



Analyticity of the ground state eigenvalues in parameters and asymptotic expansions (see Section 10 1: 

lfT5ll^[TT1l52ll63ll64ll65l. 



Existence of the ionization threshold (see Section [10]>: [50]. 
Resonance poles (see Section 10 1: [ 1 1 (see also 1 15 1). 



Self-energy and binding energy: ll78l l85l l86l l87l l56l fl9l l32l l60l IBTl l59l l58l IB1 122ll20ll. 
Electron mass renormalization: 115911791 I8"l l27ll47l 
One particle states: ||39l |40l [30l |3~T] 1471 . 
Scattering amplitudes: |[T0l . 

Semi-relativistic Hamiltonians: Il2l1 l93l l90l [84l fT05lr76ll . 
Other aspects of non-relativistic QED: 
Photo-electric effect: lfT8ll54l . 
Scattering theory: ||4T1 l42l l43l . 
Stability of matter: 
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There is an extensive literature on related models, which we do not mention here: Nelson model describing 
a particle linearly coupled to a free massless scalar field (phonons), semi-relativistic models, based on Dirac 
equation, and quantum statistics (open systems, positive temperature) models. (In the latter case, one deals 
with Liouvillians, rather than Hamiltonians, on positive temperature Hilbert spaces. The main results above 
were proved simultaneously for the QED and Nelson models and extended, at least partially, to positive 
temperatures.) 



A Hamiltonian of the Standard Model 

In this appendix we demonstrate the origin the quantum Hamiltonian H g given in ([TJ. To be specific consider 
an atom or molecule with n electrons interacting with radiation field. In this case the Hamiltonian of the 
system in our units is given by 

" 1 

H(a) = ^r( iV *i - ^ A A*i)? + <*V(x) + H f , (31) 
i=i m 

where aV(x) is the total Coulomb potential of the particle system, m is the electron bare mass, a = j^— « 

(the fine-structure constant) and A X '(y) is the original vector potential with the ultraviolet cut-off 
Rescaling x — > a~ 1 x and k — > a 2 k, we arrive at the Hamiltonian ([T}, where g :— o?l 2 and A(y) = A x (ay), 
with x(fc) := x'(a 2 fc). After that we relax the restriction on V(x) by allowing it to be a standard generalized 



n-body potential (see Subsection 2.1 1. Note that though this is not displayed, A(x) does depend on g. This 
however does not effect the analysis of the Hamiltonian H . (If anything, this makes certain parts of it simpler, 
as derivatives of A(x) bring down g.) 



In order not to deal with the problem of center-of-mass motion which is not essential in the present context, we 
assume that either some of the particles (nuclei) are infinitely heavy (a molecule in the Born-Oppenheimer 
approximation), or the system is placed in a binding, external potential field. (In the case of the Born- 
Oppenheimer molecule, the resulting V(x) also depends on the rescaled coordinates of the nuclei, but this 
does not effect our analysis except of making the complex deformation of the particle system more compli- 
cated (see [82]).) This means that the operator H p has isolated eigenvalues below its essential spectrum. The 
general case is considered below. 

In order to take into account the particle spin we change the state space the particle system to H p = 
(g)™L 2 (K 3 , C 2 ) (or the antisymmetric in identical particles subspace thereof), and the standard quantum 
Hamiltonian on U = H p ® U f , is taken to be (see e.g. l33l[34l[TU51 ) 

n l 

H ^ = S 2^ [(J ° ' ~ gA ^ x ^ + V ^ x) + H f> (32) 

3=1 

where dj := (oji, &j2, &ji are the Pauli matrices of the j— th particle and the identity operator on C 2 ™ 
is omitted in the last two terms. It is easy to show that 

[a ■ (iV a - gA x (x))} 2 - (jV x - gA x {x)f + go ■ B(x). (33) 

where B(x) := cm\A x (x) is the magnetic field. As a result the operator ([34| can be rewritten as 

- 1 

H spin = H®l + l®g22 —a 3 ■ B( Xj ). (34) 

3 = 1 m 

For the semi-relativistic Hamiltonian, the non-relativistic kinetic energy J- \p\ 2 is replaced by the relativistic 
one, ^/|p| 2 + m 2 or yj (a ■ p) 2 + to 2 . 
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B Translationally Invariant Hamiltonians 

If we do not assume that the nuclei are infinitely and there are no external forces acting on the system, 
then the Hamiltonian ([T| is translationally symmetric. This leads to conservation of the total momentum 
(a quantum version of the classical Noether theorem). Indeed, the system of particles interacting with the 
quantized electromagnetic fields is invariant under translations of the particle coordinates, x —> x + y, where 
y = (y, . . . , y) (n— tuple) and the fields, A(x) — > A(x — y), i.e. H commutes with the translations T y : 
^fe) ~~ ► e ly ' Pf ^(x + y), where Pf is the momentum operator associated to the quantized radiation field, 



P f = ^ I dkka* x {k)a x {k). 



It is straightforward to show that T y are unitary operators and that they satisfy the relations T x+y = T x T y , 
and therefore y — > T y is a unitary Abelian representation of M 3 , Finally, we observe that the group T y is 
generated by the total momentum operator, Ptot, of the electrons and the photon field: T y — e ly ' Pt °*. Here 
P to t is the selfadjoint operator on H, given by 

Ptot ■■= J2 Pi ® V + lel ® P f (35) 

i 

where, as above, pj :— —iV Xj , the momentum of the j— th electron and Pf is the field momentum given 
above. Hence [H, Ptot] = 0. 

Let T-L be the direct integral % = J R3 'HpdP 1 with the fibers Hp := L 2 (X) <g> T, where X := {x € 
R 3n | Y^i m i x i = °} - K 3 ^- 1 ), (this means that U = L 2 (R 3 , dP: L 2 (X)® F)) and define U : U e y^Hf -> 
T~L on smooth functions with compact domain by the formula 

(UV)(x',P)= [ e^ p - p ^^(x' + x cm )dy, (36) 



where x[ are the coordinates of the N particles in the center-of-mass frame and x cm — (x cm , . . . , x cm ) (n— 
tuple), with x cm = 1 m . 2~2i m i x i, the center-of-mass coordinate, so that x = x' + x cm . Then U extends 
uniquely to a unitary operator (see below). Its converse is written, for &(x', P) G L 2 (X) <E) T , as 

(U- 1 $)(x)= [ e- lx ^< p - p t1<P{x',P)dP. (37) 



The functions $(V, P) = (t/ XP) (a;', P) are called fibers of One can easiely prove the following 

Lemma B.l. The operations (f36]l and ((37]) define unitary maps L 2 (R 3n ) ® T ->• H and U ->■ L 2 (R 3n ) ® T, 
and are mutual inverses. 

Since H commutes with Ptot, it follows that it admits the fiber decomposition 

UHlf- 1 = / H(P)dP, (38) 



where the fiber operators H(P), P € M. 3 , are self-adjoint operators on T . Using a(k)e ly ' Pf =e 1J, '' Pf+fe )a(&;) 
and a*{k)e- iy - Pi = e~ ly ^ Pi ^ a* {k), we find V y e iy < p ~ p ^A x (x' + y)e lv < p - p ^ = and therefore 

A x {x)e ly -^ p - p ^ = e ly - (p - p ^A x (x - y). (39) 

Using this and ([37]), we compute H(lJ- 1 ^)(x) = J R3 e ix < p - p ^ H(P)${P)dP, where H(P) are Hamilto- 
nians on the space fibers Hp := J- given explicitly by 

H ^ = E^( P ^ f " _ e ^x(^)) 2 + Kouife') + Hi (40) 

j 1 
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where x\ = Xi — x cm is the coordinate of the i— th particle in the center-of-mass frame. Now, this hamiltonian 
can be investigated similarly to the one in ([TJ. 



C Proof of Theorem Q 

In this appendix we omit the subindex p at \ p and x p , and replace the subindex p in other operators by the 
subindex \. Moreover, we replace H — A by H. Though \ and X we deal with are projections, we often keep 
the powers x 2 and X 2 , which occur often below, having in mind showing possible generalization to x an d x 
which are 'almost (or smooth) projections' satisfying x 2 + X 2 — 1 ( see Appendix [d]). 

First we note that the relation between ip and if in Theorem 6.1 (ii) is ip = xi J , 4> — Qx(.H)fi and 
between H^ 1 and F x (iJ) _1 in (iv) is 

H- 1 = Q X (H) F X H- X Q X (H)* + xH- 1 ^ (41) 
where H x := x p Hx x and Q X (H) and Q X (H)# are the operators, given by 
Q X (H) := X - xH^Hx, 
Q#{H) := X - xHxH= l x- 



Proof of Theorem 6.1 Throughout the proof we use the notation F :— F X (H), Q := Q X (H), and 
Q# : = Q*(H). Note that (i) (0 £ p(H) «0e p(F x {H))) follows from (iv) (H- 1 exists O F x (i7) _1 
exists) and (iv) follows from ( |4T| ), so we start with the latter. 
Proof of ( pTfl ). The next two identities, 

HQ = X F and Q* H = F X , (42) 

are of key importance in the proof. They both derive from a simple computation, which we give only for the 



first equality in ( 42 1. We observe the relations 

H X = xH x +x 2 H X , and Hx = X% + X 2 H% (43) 
which follow from x 2 + X 2 — 1- Now, using the definition of the operator Q and the relations (43) , we obtain 
HQ = X H X + x 2 H X - (xH x + x 2 Hx)H- 1 xH X . (44) 
Canceling the second term on the r.h.s. with the first term in the parentheses in the third term, we see that the 



r.h.s is equal x F, which gives the first equality in (42 1 



Now, suppose first that the operator F has bounded invertible and define 

R := QF- l Q* + xH x x x- (45) 

Using ( |42) and ( |43] l, we obtain 

HR - HQF~ X Q* + (xH x + x 2 Hx)H- 1 x (46) 



XQ* + X 2 + X 2 HxH= x x 



= x 2 + x 2 = 1, 

and, similarly, RH = 1. Thus R = H~ l , and ( |4T) holds true. 

Conversely, suppose that H is bounded invertible. Then, using the definition of F and the relation x 2 + 
X 2 = 1, we obtain 

F X H- X x = XH X 2 H~ x x- xHxH^xHx 2 H~ X X (47) 
= X Hx 2 H-'x- X HxH^xHH- 1 x+ xHx H^xHt B~ X X 



X Hx 2 H~ l x+ xHx 2 H-\ = X 2 - 
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Similarly, one checks that \ H 1 \F = 1. Thus F is invertible on Ran \ with inverse F 1 = x H 1 x- 
Proof of (ii) (Hip = XtP F p (H - A) tp = 0). Ifip eH\ {0} solves i?^ = then (p2ll implies that 



Furthermore, by (|43]l, = x H tp = H x x*P 



'X"0 = X ^ 



FxV> = Q # i?^ = 0. 

xHx 2 i/j, and hence 



xH^xHx 2 ^ = X 2 ^ 



Therefore, V> 7^ implies x4> 7^ 0. 

If tp G Ran x \ {0} solves Ftp = then the definition of Q implies that 

xQtp = xf = f, 



(48) 



(49) 



(50) 



which implies that Qtp 7^ provided tp 7^ 0. 

Proof of (iii) (dim Null(tf-A) = dimNullF p (??-A)). By(i), dimNulli? = is equivalent to dim Null/ 1 = 
0, assuming that H £ D(F). We may therefore assume that Nulli? 7^ and NulLF 7^ are both nontrivial. 



Eq. (49 1 shows that x ■ Nulli? — > Null/ 1 is injective, hence dimNulli/ < dim NulLF, and Eq. ( 50 » shows 



that Q : Nulli 7, — > Nulli? is injective, hence dimNullii > dim NulLF. This establishes (iv) and moreover 
that x ■ Nulli? -> NulLF and Q : NullF -> Null/? are actually bijections. □ 



D Smooth Feshbach-Schur Map 
D.l Definition and isospectrality 

We define the smooth Feshbach-Schur map and formulate its important isospectral property Let x, A 7 be a 
partition of unity on a separable Hilbert space ??, i.e. x and X are positive operators on ?? whose norms are 
bounded by one, < x, x < 1, and x 2 + X 2 = 1- We assume that x and x are nonzero. Let r be a (linear) 
projection acting on closed operators on 7? with the property that operators in its image commute with x an 
r(l) = 1. Let r := 1 — t and define 

H T<X# :=r(H) + x*r(H) X * , (51) 

where x^ stands for either x or % 

Given x and r as above, we denote by D T:X the space of closed operators, H, on ?? which belong to the 
domain of r and satisfy the following three conditions: 

(i) t and x (and therefore also f and x) leave the domain D(H) of i? invariant: 

D{t{H)) = D(H) and X D{H) C D(H), (52) 

(ii) 

i?T,x i s (bounded) invertible on Ranx, (53) 

(iii) 

r(i?)x and xt(H) extend to bounded operators on ??. (54) 

(For more general conditions see lP7l l5"Tll .) 

The smooth Feshbach-Schur map (SFM) maps operators from D TtX into operators on 7? by 

F T)X (i?) := H + x.W X - X^xH^xWx, (55) 

where Ho := t(H) and W := t(H). Note that Ho and W are closed operators on ?? with coinciding 
domains, D(H ) = D(W) = D(H), and i? = i? + W. We remark that the domains of xWx, xW% Ht, x < 
and i? T ,x a U contain D(H). 

Define operators Q T:X (H) := x ~ xH~^xWx and Q* X {H) := X ~ xWxH~^x- The following 
result ([7]) generalizes Theorem[6T]above; its proof is similar to the one of that theorem: 
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Theorem D.l (Isospectrality of SFM). Let < x < 1 and H G D Ta be an operator on a separable Hilbert 
space H. Then we have the following results: 

(i) H is bounded invertible on W if and only if F T , X (H) is bounded invertible on Ranx- In this case 

H- 1 = Q T>X {H) F TiX {H)~ l Q T , X {H)* + xH^ l % (56) 
F^H)- 1 = X H~ X X + XriH)- 1 ^. (57) 

(ii) If i/j £ FL \ {0} solves Hijj = then <p := xi> € Ranx \ {0} solves F T<X (H) ip = 0. 

(iii) If cp e Ranx \ {0} solves F T>X (H) p = then ip := Q TyX (H)(p eH \ {0} solves = 0. 

(iv) The multiplicity of the spectral value {0} is conserved in the sense that dim Null H = dim Null F TyX (H). 
We also mention the following useful property of F T , X : 

H is self-adjoint =>■ F T X (H) is self-adjoint. (58) 

D.2 Transfer of local decay 

We have shown above that the smooth Feshbach-Schur map is isospectral. In fact, under certain additional 
conditions it preserves (or transfers) much stronger spectral property - the limiting absorption principle (LAP) 
([46 1), which is defined as follows. Let A C i be an interval, v > and B, a self-adjoint operator. We say 
that a C 1 family of self-adjoint operators, s.t. H (A) € D T X has the (A, B, 9) limiting absorption principle 
(LAP) property iff 

Um (B)- B (H(X) - iey 1 {B)- e exists and e C 1/ (A). (59) 

Usually LAP holds for v < 6 - \. One can show that the LAP implies the local decay property (see, e.g. 
l99l . vol III; recall that the definition of the local decay property is given in SectionfTO]). 

Theorem D.2. Let A C K and , VA e A, H(X) be a C 1 family of self-adjoint operators, s.t. H(X) e D r . x . 
Assume that there is a self-adjoint operator B s.t. 

ad-^(A) is bounded and differentiable in A, Vj < 2, (60) 

where A is one of the operators x, Xi Xr(H(X)), f(H(X))x, d^(xH TtX (X)~ 1 x) , k = 0,1. Then, for any 
< v < 1 and < 6 < 1 and in the operator norm, we have 



Urn (B)- e (F T . X (H(X)) - ie)^ 1 {B)-° exists and e C U (A) (61) 
lim + (B)- e (H(X) -iey 1 (B)- e exists and e ^(A). (62) 

This allows one to reduce the proof of the LAP for the original operator, H — A, to the proof of this 
property for a much simpler one, lZp(H — A). 
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D.3 Transfer of analyticity 

Theorem D.3. Let A be an open set in C and H(X), A G A, a family of operators with a fixed domain, which 
belong to the domain of F TX . Assume H(X) and t(H(X)), with the same domain, are analytic in the sense 
of Kato (see e.g. [99 1, vol IV). Then we have that 

• F TtX (H(X)) is an analytic in A G A family of operators. 

Proof. Since that H(X), H (X) := t{H{X)) and W{X) X and xW{X), where W(X) := t(H(X)), are 
analytic in A G A, we see from the definition of the smooth Feshbach-Schur map F TX in ( |55] > that F TX (H(X)) 
is analytic in A G A, provided x p H (A)~i x P is analytic in A. The analyticity of the latter family follows by 
the Neumann series argument. □ 

One can generalize the above result to A's which are open sets in a complex Banach space. Recall that a 
complex vector-function / in an open set A in a complex Banach space VV is said to be analytic iff it is locally 
bounded and Gateaux-differentiable. One can show that / is analytic iff G W, f(H + t£) is analytic in 
the complex variable t for |r| sufficiently small (see |23, 66 1). Furthermore if / is analytic in A and g is an 
analytic vector-function from an open set ft in C into A, then the composite function / o g is analytic on f2. 



E Mass Renormalization 

As the free electron is surrounded by virtual 'soft' photons its effective (inertialO mass is greater than the 
value ('bare' mass) entering its Hamiltonian. One calls this electron mass renormalization. We begin with 
analyzing the definition of (inertial) mass in Classical Mechanics. Consider a classical particle with the 
Hamiltonian h(x, k) :— K(k) + V(x), where K(k) is some function describing the kinetic energy of the 
particle. To find the particle mass in this case we have to determine the relation between the force and 
acceleration at very low velocities. The Hamilton equations give x = d^K and k = F, where F = —d x V 
is the force acting on the particle. Assuming that K has a minimum at k = and expanding dkK(k) around 
0, differentiating the resulting relation x = K"(0)k, where K"(0) is the hessian of K at k = 0, w.r. to 
time and using the second Hamilton equation, we obtain x. = K" (0)F(x). This suggests to define the mass 
of the particle as m — AT"(0) _1 , i.e. as the inverse of the Hessian of the energy, in the absence of external 
forces, as a function of of momentum at 0. (K(k) is called the dispersion relation.) We adopt this as a general 
definition: the (effective) mass of a particle interacting with fields is the inverse of the Hessian of the energy 
of the total system as a function of of the total momentum at 0. 

Now, we consider a single non-relativistic electron coupled to quantized electromagnetic field. Recall that 
the charge of electron is denoted by — e and its bare mass in our units is m. The corresponding Hamiltonian 
is 

H ■= ^—(iV x ®lf - eA x (x)) 2 + latoHf, (63) 
2m 

acting on the space L 2 (M 3 ) <Ei T = H pai t < ^>'H-f- It is the generator for the dynamics of a single non-relativistic 
electron, and of the electromagnetic radiation field, which interact via minimal coupling. Here recall A x (x) 
and Hf are the quantized electromagnetic vector potential with ultraviolet cutoff and the field Hamiltonian 
and are defined in |2]) and ([3]l 

The system considered is translationally invariant in the sense that H commutes with the translations, 

Ty, 

TyH = HTy, 

which in the present case take the form 

T y : -> e vy Pt ^{x + y), (64) 
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This as before leads to H commuting with the total momentum operator, 

Ptot := Pel® 1/ + l e i ® Pf, (65) 

of the electron and the photon field: [H, Ptot) — 0. Here P c \ := — iV x and Pf = J^x J dkk a* x (k)a\(k) are 
electron and field momenta. Again as in Appendix [B] this leads to the fiber decomposition 

UHlf- 1 = / H(P)dP, (66) 



where the fiber operators H(P), P € M 3 , are self-adjoint operators on T. Using a(k)e lx ' Pi =e lx '( Pl+kS, a(k) 
and a*{k)e' lx - Pt = e~ lx < Pi ~ k) a*(k), we find V x e lx < p ~ Pi) A x {x)e lx < p - Pf) = and therefore 

A x (x)e lx < p - Pf ) = e lx < p - Pt )A x (0). (67) 

Using this and (J37j>, we compute i/(J7" 1 $)(a:) = J R3 e ix < p - p ^ H(P)$(P)dP, where H{P) are Hamilto- 
nians on the fibers Hp := J- given explicitly by 

H(P) = J- (P - P f - e A x f + Hi (68) 



2m y ~ 

where A x := A X (Q). Explicitly, A x is given by 



A x - E / dfc S { o A (fc) + al(k) }. (69) 



Consider the infimum E(P) := inf a(H(P)) of the spectrum of the fiber Hamiltonian H(P). Note that 



for e = 0, £(P)| e=0 =: £fc(P) is the ground state energy of if (P) := H{P)\ e=0 = ^(P- Pf) 2 + fff 

ip|2 

with the ground state ^1 and is E (P) = The renormalized electron mass is defined as 

E{P) = + 0(\P\ 3 ) 

where the left hand side is computed perturbatively up the second order in the coupling constant (charge). 
Provided that E(P) is spherically symmetric and C 2 at P = 0, and therefore, in particular, d\p\E(Q) = 0, 
we define the renormalized electron mass at zero total momentum as 



1 



a 2 pl p(o) 



The kinematic meaning of this expression is as follows. The ground state energy E{P) can be considered as 
an effective Hamiltonian of the electron in the ground state. (The propagator exp(— itE(P)) determines the 
propagation properties of a wave packet formed of dressed one-particle states with a wave function supported 
near p = - which exist as long as there is an infrared regularization.) The first Hamilton equation gives the 
expression for the electron velocity as 

v = dpE(P) . 

Expanding the right hand side in P we find v = Hess E(0)P + 0(P 2 ), where 

(Hess^(P)) .. = (5 t] - m + ^df Pl E(P) (70) 
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is the Hessian of E(P) at P £ M 3 (given that E(P) is spherically symmetric, and C 2 in \P\ near P — 0). It 
follows from (PTol and the fact d\ P \E(Q) = that Hess£(0) = d? p ,E{Q) 1, so that 

v = df Pl E(0,a)P + 0{P 2 ) . 

This suggests taking (^^^(O))" 1 as the renormalized electron mass at P = 0. 
The following result is proven in |8, 27, 30 1: 

Theorem E.l. For any P, s.t. \P\ < |, the infimum of the spectrum E(P) — inf spec(-ff (P)) is twice 

< 1 + eg 2 for some c > 0. 



A presentation of the leading order calculations can be found in [79|. 

Remark 1. The estimate 1 < m ren < 1 + eg 2 reflects the fact that the mass of the electron is increased 
by interactions with the photon field. 



F One-particle States 

First we note that for e = 0, the one-particle states of the Hamiltonian Hq := H\ e -o = \P c i\ 2 + H{ are the 
generalized eigenfunctions 

e- iP x ® <!. (71) 

|p|2 

corresponding to the spectral points E(P) = -^[- This corresponds to the true ground state ft of the fiber 
Hamiltonians Hq(P) := H(P)\ e=0 = oa-P — Pf) 2 + Hi- The generalization of such a state for the 
interacting system would be the ground state of H(P), if it existed. However, we have 

Theorem 1. H(P) has a ground state if and only if P = 0. 



To define one particle states for the interacting model, we first introduce IR regularization, K a £ ( [0, n] ; R+ ) 
is assumed to be a smooth cutoff function obeying lim^^o K " 1 = 1 • The corresponding Hamiltonian is 

1 2 
H a := — ( -iV x ® If + eA a (x)) + l el ®H f (72) 

with the quantized electromagnetic vector potential subjected, besides the ultra-violet cutoff, also infrared 
regularization, 

M*) = T,J dk ^07Y { ^ ® a ^ + h - c - >• (73) 
Let 5 :={PeR 3 \ \P\ < 1/3}. 

Theorem 2. For P e S and for any a > 0, the infimum of the spectrum E a (P) = inf spcc(iJ CT (P)) is a 
simple eigenvalue. 



Remark 2. The upper bound on \P\ of ^ is not optimal, but we note that, For E(P) to be an eigenvalue, 
\P\ cannot exceed a critical value P c < 1 (corresponding to the speed of light). As \P\ — > P c , it is expected 
that the eigenvalue at E{P) dissolves in the continuous spectrum, while a resonance appears. This is a 
manifestation of a phenomenon analogous to Cherenkov radiation. 
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Let ^o-(P) e T denote the associated normalized fiber ground state, ||^ CT (P)||jr = 1, for P E S, 

H„(P)* a (P) = E„{P) *„(P) . (74) 

The vector &(P, <r) is an infraparticle state, describing a compound particle comprising the electron together 
with a cloud of low-energy (soft) photons whose expected number diverges as a — > 0, unless p = 0. For 
f £ 5, we introduce the Weyl operators 

^ (P) W:=e D W- D *W (75) 

where 

: = E / dfcG A (fc,p) e - 4fea; aA (fc), (76) 

with 

G A (M := Ml*l) |fc|1/2 (| fc |_ V ^( P ).fe) ' (7?) 

(Here and in the sequel, we will use the abbreviated notation VE a (p) = V p E(p) .) We observe that they 
commute with the total momentum operator, 

[Ptot, W VEAp) (x)} = 0. (78) 

To see this, we note that [P tot , D(x)] = = [P to t, D* {x)\. Indeed, we have that 

Ptot Y, I dkG x {k,p)e- lkx a x {k)i;{x) 
x J 

= I dkG x (k,p) (-iV x + P f ) e~ ikx a x {k)i>{x) 

x J 

= Y [ dkG x (k,p)e- tkx a x (k)(-i\7 x + k + P f -k)iP(x) 
x 

= [ dkG x (k,p)e- ikx a x (k)P tot ^(x). (79) 

Accordingly, we infer that W VEAp -j(x) — exp[D(x) — D*(x)} commutes with Ptot- 
Furthermore, we observe that 

W VEAp) {x)e^~ p ^ = e i(p-Pf)*W VEAp) (80) 

holds. Here and in what follows, we will use the abbreviated notation 

W VEAp) = W VEAp) (x = 0). (81) 

We define the maps 

:= J d P W VEAp) (x)e^- p ^j(p) 

= J dpe^- p ^W VEAp) d(p)- (82) 

Likewise, 

:= j ' dpW* WEAp) {x)e^- p ^m. (83) 
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The associated Bogoliubov-transformed Hamiltonian is given by 

K„ := (WH a W*). (84) 
We also introduce the Bogoliubov-transformed fiber Hamiltonians 

K a (p) := Wve Ap ) H a (p) W^ EAp) . (85) 

Then,we observe that 

K a = (yvH a yv*){x) 

W VEAp) (x) e^- p ^H(p)e-^- p ^ dP Ptot (p) W^ EAp) (x) 
W VEAp) (x) e ^- p f>H(p)e-^- p f> W* E<r(p) (x) dP Ptot { P ) 
e^ p - p f>K x {p)e-^ p - p ^ x dPp tot (p). (86) 



In particular, we have that 

W(I^)= W), (87) 

as can be readily verified. Defining 

*„(p) := W VM p)*«r(p), (88) 

we obtain 

K a (p)^ a (p) = E a (p)<5> a (p). (89) 
The following result is proven in (30). 

Theorem 3. For any P G S, the ground state eigenvector $ CT (P) of K a (P) converges strongly in T: 
c(>(p) : = lim^-^o ®<t(P) exists in T. 



G Pull-through formulae 

In this appendix we prove the very useful "pull-through" formulae (see 0~3]) 



and 



a(k)f(H f ) = f(H f + u{k))a{k) 



f(H f )a*(k) = a*(k)f(H f +uj(k)), 



(90) 



(91) 



valid for any piecewise continuous, bounded function, /, on K 
a{k), a*{k), one proves relations (90l- (91 1 for f(H) = (H f 



First, using the commutation relations for 



Weierstrass theorem, one can extend (90 1- (91 1 from functions of the form /(A) = (A — z) 
the class of functions mentioned above. 



G C/M+. Then using the Stone- 



Z G 



to 
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H Supplement: Creation and Annihilation Operators 

Let \] be either L 2 (R 3 , C, d 3 k) or L 2 (R 3 ,C 2 , d 3 k). In the first case we consider f) to be the Hilbert space of 
one-particle states of a scalar boson or phonon, and in the second case, of a photon. The variable k <E M 3 is 
the wave vector or momentum of the particle. (Recall that throughout these lectures, the propagation speed 
c, of photon or photons and Planck's constant, H, are set equal to 1 .) The Bosonic Fock space, T, over f) is 
defined by 

oo 

J:=05„(,®", (92) 

n=0 

where S n is the orthogonal projection onto the subspace of totally symmetric n-particle wave functions 
contained in the n-fold tensor product h®" of h; and 6>of)®° : = C. The vector £1 := (1, 0, ...) is called the 
vacuum vector in JF. Vectors $ £ J can be identified with sequences (ip n )^Lo of n-particle wave functions, 
which are totally symmetric in their n arguments, and ipo g C. In the first case these functions are of the 
form, ip n (ki, . . . , k n ), while in the second case, of the form ip n (ki, X%, . . . , k n , X n ), where Xj G {—1, 1} are 
the polarization variables. 

In what follows we present some key definitions in the first case, limiting ourselves to remarks at the end 
of this appendix on how these definitions have to be modified for the second case. The scalar product of two 
vectors ^ and <£> is given by 

oo „ n 

k n ) <Pn(kl, . . . ,kn)- (93) 

n=0 ^ j=l 

Given a one particle dispersion relation uj(k), the energy of a configuration of n non-interacting field 
particles with wave vectors k%, . . . , k n is given by Y^j=i We define the free -fie Id Hamiltonian, Hf, 

giving the field dynamics, by 

n 

{H f V) n (h,...,k n ) = ($^w(%))Vn(*l,...,fc*), (94) 

j=l 

for n > 1 and (Hf^>) n = for n = 0. Here ^ = {ipn)n°=o ( t0 De sure mat tne r -h- s - makes sense we 
can assume that ip n = 0, except for finitely many n, for which ip n (kx, . . . ,k n ) decrease rapidly at infinity). 
Clearly that the operator Hf has the single eigenvalue with the eigenvector £1 and the rest of the spectrum 
absolutely continuous. 

With each function tp € L 2 (M. 3 , C, d 3 k) one associates an annihilation operator a((p) defined as follows. 
For ^ = (V-'n)^Lo "= -F w i tn tne property that ip n = 0, for all but finitely many n, the vector a(ip)^f is defined 
by 



(a((p)^) n (k u ...,k n ) := Vn + 1 / d 3 k cp[h) ip n +i(k, fei, . . . , k n ). (95) 



These equations define a closable operator a((p) whose closure is also denoted by a{<p). Eqn (??) implies the 
relation 

a(<p)n = 0, (96) 

The creation operator a*(ip) is defined to be the adjoint of a(ip) with respect to the scalar product defined in 
Eq. (|93|). Since a(ip) is anti-linear, and a* (tp) is linear in ip, we write formally 



a(ip) = d d k(p(k)a(k), a*(ip) = d 3 k <p(k) a*(k), (97) 



where a(k) and a*(k) are unbounded, operator-valued distributions. The latter are well-known to obey the 
canonical commutation relations (CCR): 

[a # (fc), a*(k')] = 0, [a(jfc), a* (A;')] = S 3 (k-k'), (98) 
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where a # = a or a* . 

Now, using this one can rewrite the quantum Hamiltonian Hf in terms of the creation and annihilation 
operators, a and a*, as 

H f = J d 3 k a* (ft) u(k) o(jfe), (99) 

acting on the Fock space J 7 . 

More generally, for any operator, t, on the one-particle space L 2 (M. 3 , C, d 3 fc) we define the operator T on 
the Fock space T by the following formal expression T := fa* (k)ta(k)dk, where the operator t acts on the 
A:— variable (T is the second quantization of t). The precise meaning of the latter expression can obtained by 
using a basis {4>j} in the space L 2 (M 3 , C, d 3 k) to rewrite it as T := 2~2j J a* ((j>j)a(t*<f>j)dk. 

To modify the above definitions to the case of photons, one replaces the variable k by the pair (ft, A) and 
adds to the integrals in k sums over A. In particular, the creation- and annihilation operators have now two 
variables: a* (ft) = a#(fc, A); they satisfy the commutation relations 

[a*(k),a*(k')] = 0, [a x (k), a* x ,(k')] = 6 x , y 6 3 (k ~ k r ). (100) 
One can introduce the operator- valued transverse vector fields by 

a#(ft):= J2 ex(k)a*(k), 
ag{-i,i} 

where e^(fc) = e(fc, A) are polarization vectors, i.e., orthonormal vectors in M 3 satisfying ft • e^(fc) = 0. 
Then, in order to reinterpret the expressions in this paper for photons, one either adds the variable A, as was 
mentioned above, or replaces, in appropriate places, the usual product of scalar functions or scalar functions 
and scalar operators by the dot product of vector-functions or vector-functions and operator-valued vector- 
functions. 
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